Random global coupling induces synchronization and nontrivial collective behavior in 

networks of chaotic maps 
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The phenomena of synchronization and nontrivial collective behavior are studied in a model of 
coupled chaotic maps with random global coupling. The mean field of the system is coupled to 
a fraction of elements randomly chosen at any given time. It is shown that the reinjection of the 
mean field to a fraction of randomly selected elements can induce synchronization and nontrivial 
collective behavior in the system. The regions where these collective states emerge on the space of 
parameters of the system are calculated. 

PACS numbers: 05.45.-a, 89.75.Kd 



There is much current interest in the investigation of 
collective properties of complex networks of interacting 
nonlinear elements As a model to study some mini- 
mal conditions for the emergence of collective behavior in 
a chaotic network, we consider the coupled map system 



b t+i 
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with probability p, 
with probability 1 — p, 
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(i = 1,2, ... ,N; N = system size) gives the 
state of the ith element at discrete time t; e is the cou- 
pling strength, / is a map defining the local dynamics, 
and 
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is the instantaneous mean field of the system that pro- 
vides a global coupling. The parameter p is the proba- 
bility of connection of an element to the mean field at 
time t. Thus the average fraction of globally connected 
elements at any given time is p. 

As local dynamics in Eq. |T]) we shall consider the 
logarithmic map f(x) — b + log \x\ Q, where b is a real 
parameter and x 6 (—00,00). This map exhibits robust 
chaos, with no periodic windows and no separated chaotic 
bands, in the parameter interval b G [— 1, 1]. 

A synchronized state at time t is defined by the con- 
dition x\ = xj,yi,j, in which case the dynamics is de- 
scribed by the single map Xt+i = f(xt)- The synchro- 
nization of the elements in the network can be character- 
ized by the time-average (a) of the instantaneous stan- 
dard deviations at of the distribution of site variables x\ , 
denned as 
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Figure 1 shows the quantity (<j) (right vertical axis) as 
a function of the probability p, for fixed values of b and 



e. There is a threshold value p c ~ 0.75 at which (a) = 
(within a precision of 10 -8 in our calculations), indicat- 
ing that the elements are synchronized. The range of 
p where chaotic synchronization takes place is indicated 
by the label S. For p = 1, the system Eq. {1} becomes 
a globally coupled map network [3j which is known to 
synchronize. However, Fig. 1 reveals that the reinjection 
of a global coupling function to a fraction of randomly 
selected elements in the system is enough to achieve syn- 
chronization. The critical value p c for the emergence of 
synchronization depends on the parameters of the sys- 
tem, as shown in Fig. 2. 




FIG. 1: Left vertical axis: bifurcation diagram of ht as a func- 
tion of p. For each value of p, the mean field was calculated 
at each time step during a run starting from random initial 
conditions on the local maps, uniformly distributed on the 
interval [—8,8], after discarding the transients. The regions 
where collective states occur are labeled T (turbulent), NTCB 
(nontrivial collective behavior), CB (chaotic bands), S (syn- 
chronization). Right vertical axis: (a) vs p, continuous line. 
Fixed b = -0.7, e = 0.4; size N = 10 5 . 
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The instantaneous mean field of the system h t can 
characterize more complex collective behaviors arising in 
the system Eq. (JT]). When b 6 [—1,1], the elements in 
the network are chaotic and desynchronized. However, 
for some parameter values, ht reveals the existence of 
global periodic attractors. Figure 1 shows the bifurca- 
tion diagram of h t (left vertical axis) as a function of p, 
for fixed b and e. In this representation, collective peri- 
odic states at a given value of the parameter p appear 
as sets of small vertical segments which correspond to 
intrinsic fluctuations of the periodic orbits of the mean 
field. 

In the region labelled by T (turbulent) in Fig. 1, h t fol- 
lows a normal statistical behavior around a mean value 
(a collective fixed point) with fluctuations reflecting the 
averaging of N completely desynchronized chaotic ele- 
ments. Increasing the probability of connection p induces 
a transition to periodic collective states occurring in the 
chaotic range of the local dynamics: a pitchfork bifurca- 
tion takes place from a collective fixed point to a collec- 
tive period-two state (a state for which the time series of 
h t alternatingly moves between the corresponding neigh- 
borhoods of two separate values). Collective states of 
higher periodicity arise by further increasing p: global 
periodic attractors of period 2, 4, 8, and 16 are possible 
in this system. The amplitude of the collective periodic 
motions manifested in h t do not decrease with an in- 
crease in the system size N. As a consequence, the vari- 
ance of ht itself does not decay as TV -1 with increasing 
N but rather it saturates at some constant value related 
to the amplitude of the collective period. This is a phe- 
nomenon of nontrivial collective behavior, where macro- 
scopic quantities in a spatiotemporal dynamical system 
exhibit ordered evolution coexisting with local chaos [ij. 
Note that the emergence of collective periodic behavior in 
this system cannot be attributed to the presence of peri- 
odic windows in the local dynamics since the logarithmic 
map possesses robust chaos for b <E [—1,1]. Figure 1 indi- 
cates with the label NTCB the interval of p where non- 
trivial collective behavior arises in this system. Before 
crossing the boundary of the synchronization region, the 
collective states described by h t take the form of chaotic 
bands. These states are labelled CB (collective bands) 
and they consist of the motion of chaotic elements that 
maintain some coherence. 

Figure 2 shows the regions where the different collec- 
tive states of the system Eq. ([J) occur on the space of 



parameters (p, e) . These regions are separated by stabil- 
ity boundaries which have been numerically calculated. 
For p = 1, Fig. 2 yields the intervals of stability of the 
collective states S, NTCB and CB corresponding to glob- 
ally coupled logarithmic maps and which have been pre- 
viously found [a]. Figure 2 shows that those collective 
states can also emerge when a fraction of randomly se- 
lected maps are connected to the mean field and appro- 
priate values of the coupling strength are employed. 
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FIG. 2: Boundaries on the parameter plane (p, e) separating 
different collective states of the system. Labels correspond to 
those in Fig. 1. Fixed b = -0.7. Size N = 10 5 . 



In summary, we have shown that a reinjection of a 
global coupling function to a fraction of randomly se- 
lected chaotic elements in a dynamical network can in- 
duce synchronization and nontrivial collective behavior 
in the system. This procedure may have practical appli- 
cations in the control of spatiotemporal systems. These 
results may be relevant in some biological and social con- 
texts where the global information is often available only 
to a portion of agents in those systems. 
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